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of the t-channei partiai waves, and corresponding eigenfunctions of the BFKL equation in ieading iog(l/x) 
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1. Introduction 


In our paper [1] we solved the BFKL equation with a massive gluon in the framework of the Higgs model 
numerically. Such an equation arises in the electroweak theory with zero Weinberg angle ( see Ref. [2]). 
However, for us, this gauge invariant theory provides an instructive example of a model in which the 
scattering amplitude has correct large impact parameter behavior ( scattering amplitude oc exp {—mb) at 
large 5), but we still have the unitarity problem as the scattering amplitude increases as at high energies. 
Therefore, this model is a perfect training ground to study how the correct b behavior, can influence the 
resolution of the unitarity problem in the framework of the CGC/saturation approach [3-6]. This approach 
leads to a partial amplitude smaller than unity, as required by unitarity constraints. However, it generates 
a radius of interaction that increases as power of energy [7-10], leading to the violation of the Froissart 
bound [11]. 

Another facet of this model, is that it is a possible candidate for an effective theory equivalent to 
perturbative QCD, in the region of distances (r) shorter than 1/m, where m denotes the gluon mass. 
Indeed, for r <C 1/m, the fact that gluon has a mass, is not essential while at r ~ 1/m similar correlation 
functions arise from fixing or eliminating Gribov’s copies [12] (see Refs. [13-15]). We wish to stress that 
a gauge theory with the Higgs mechanism, leads to a good description of the gluon propagator, calculated 
in lattice approach [16] with m = 0.54 GeR. Therefore, a plausible scenario is that the Higgs gauge theory 
describes QCD in the kinematic region r < Ijm, while for r ~ I/Aqqd > 1/m the non-perturbative QCD 
approach takes over, and leads to the confinement of quarks and gluons, which is missing in the theory 
with a massive gluon. 

We found [1] that the spectrum of the massive BFKL equation in w-space for t = 0, coincides with 
that of massless one [17,18]. The simple parametrization of the eigenfunctions have also been obtained in 
Ref. [1]. In this paper we propose using the solution to the BFKL equation with massive gluons which has an 
advantage of being simple and semi-analytic. Having this solution in hand, we are able to progress to more 
difficult problems, e.g. a generalization of the main equations of the CGC/saturation approach [6,19,20]. 

The equation for the scattering amplitude in leading ln(l/x) approximation of perturbative theory, 
has been on the market for some time [17], and is schematically shown in the Fig. Its kernel is of the 
form [1,17] 


K {qi,q2Wi,q2) 
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( 1 . 1 ) 


For q = 0, the kernel Q simplihes considerably, and yields a homogeneous BFKL equation for the 
Yang-Mills theory with the Higgs mechanism. 
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Figure 1: The massive BFKL equation (Fig. |^a) and its kernel (Fig. ^b) 


where qi = q2 = p, q'l = q2 = p' and co{p) is the gluon Regge trajectory given by 
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Examining the rotationally symmetric solution, the kernel can be integrated over the azimuthal angle </>. 
Introducing the new variables 


K = 


p^ 




/2 


E = -^- d5 = (1.4) 

as vr 

and changing the notation of the wave function f{p) to we obtain the one-dimensional BFKL 

equation 
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where the kinetic energy is given by 


dn' (f)E {k') 


N^ +1 1 

0 a/(«: - Ac')^ + 2{k -b ac') + 1 2iV2 i 


+ 


f 


(f)E {k') dn' 

k’ + 1 


T(k) = 


K + 1 \/ K + 4 + 

In 


,(1.5) 


( 1 . 6 ) 


y/K\^ K + 4 y/K + 4: — y/R 
In this paper we will deal mostly with the wave function in Y representation defined by a Fourier transform, 

r-6+ioo 


4>{Y,k) = J 

For 4/ (Y, k) , the Eq. (|1.5D takes the form 
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For completeness of presentation, we recall that the massless BFKL equations have the following form 
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The way of regularization at e —?> 0 stems directly from Eq. (O) considering small masses em^ instead of 
Eq. ( |1.9|) can then , be re-written with a different way of regularization, for example 

(1.9) 

e^U \^/ Jo — K|“re 

where 0 (z) is the unit step function which is equal 1 for z > 0 and 0 for z < 0. At large k solutions of 
Eq. ( ll.SD , should coincide with the eigenfunctions of Eq. (O) which are well known [17,18] and have the 
form 


^ ~ ^ with E( 7 ) = -x{l) = V’(7) + V’(l- 7 ) - 2V^(1), (1-10) 


(j)E {n) - 

where V'(x) = r'(x)/r(x) is a digamma function. Two eigenfunctions, (k) (x and 4>e («)) on 
describe the states with the same energy. For 7 = ^ -|- w these eigenfunctions are normalized and form a 
complete set of functions. 

2. Semi-classical approach: generalities and eqnations 
2.1 The main qualitative features of the solution 


For completeness of presentation, we start discussing the solutions to Eq. (|1.5|) , repeating the key qualitative 
and general features of solutions that have been discussed in Ref. [1]. The first one has been mentioned in 
the previous section: at large values of k, the eigenfunctions 4>e (k) should approach the eigenfunctions of 


the massless BFKL equations (1.10). 

The behavior of the solutions at small values of k, is easier to understand by rewriting Eq. (1.5) in the 
coordinate representation. Using 
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where Jq (z) and Kq {z) are the Bessel and MacDonald functions [21], we can re-write Eq. (1.5) in the 
form 
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In (2.3) we introduced a shorthand notation P for the projector onto the state ~ mr j{p^ + mr) 

f d'^p' 4>{p') in coordinate representation n < \ f d^p' 4>{p') 
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The behavior of Eq. (1.10) at large p, translates into the following behavior at short distances 


„ , . r<l/m 
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To understand the behavior of the solutions at large distances, we should distinguish the two cases, when 
the wave function’s decrease is slower than , and when the decreases are faster In the former 


case, as we may see from Eq. (2.3,2.4), we may neglect the contact term and term ~ V{r), and Eq. 
degenerates into 


TfE = E{'y)i^) — 
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In 
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where k = — V^. 

The eigenfunctions of Eq. (|2.6| ) have a form 

fE = E{u)i^) ~ for a < 0; and (r) ~ for a > 0; (2.7) 


From Eq. (^) we can see that parameters a and 7 are correlated, viz. 

E = T{-a) = -x(7) (2.8) 

The solution to this equation is shown in Fig. The Fourier image of l/(p^ + in the coordinate 



representation is 
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{k + a)^“'>' 
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(2.9) 


The function 1 /(k + a)^ describes both the short distance (|2.5|) and long distance (^.71) behaviors. In [1] 
we demonstrated that the ground state with the minimal energy is reached with this class of functions. 
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2.2 Semi-classical approach: generalities 
2.2.1 The method of steepest descent 


For massless BFKL, a general solution has a form 

re+ioc 




( 2 . 10 ) 


where 4>in ( 7 ) should be found from the initial conditions at T = 0. For the massive case, we will look for 
solutions of Eq. (1.8) of the analogous form 


iY,l) = r‘“ ( 7 ) 

J e—ioQ 




( 2 . 11 ) 


where we introduced a new variable I = In (k + a), which is stable in the small-K limit, and 4>in ( 7 ) is 
fixed by the initial conditions at T = 0. In Eq. ( ^.ll| ) we can take the integral over 7 using the method 
of steepest descent, which is equivalent to a search of semiclassical solution to Eq. dU). The equation for 
the saddle point takes the general form 


dE (7SP iY,l),l) 
dj 


Y + I = 0. 


( 2 . 12 ) 


If 4 >in (7) in the integrand of (2.11) is a smooth function, we can replace it with its value at 7 = 7sp- The 
integral over 7 in the vicinity of a saddle-point (2.12|) yields 
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The omission of higher-doer corrections in Eq. (2.13|) is justified due to smallness of the parameter 
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The smoothness of the initial function (pin ( 7 ) implies a condition 
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1 d^E{-ysp{Y,l),l) 
Tp 
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^ l7=7SP 


(2.15) 


and determines the kinematic region of applicability of the semiclassical approximation. As we will demon¬ 
strate below, both conditions are satisfied for sufficiently large y S> 1. 
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2.2.2 Solution with method of characteristics 


The method of characteristics for a partial differential equation (PDE) corresponds to a reduction of the 


PDEs 


F(xi, . . .,Xn,U,pi, . . . ,Pn) = 0, 



(2.16) 


to a system of ordinary differential equations (ODE) for characteristic lines along which the PDE converts 
into an ordinary differential equation. These characteristics satisfy the Lagrange-Charpit equations 


Xj _ Pi _ u 

^Pi ^Xi T Xy Pi^Pi 


(2.17) 


A n instructive example familiar from a classical mechanics, is a Hamilton-Jacobi equation, for which 
the characteristic lines correspond to trajectories of particles which are solutions of the Newtonian equations 
of motion. A detailed discussion of the method is beyond the scope of the present paper and can be found 
in the literature (see e.g. textbooks [22,23]). 

A direct application of the method of characteristics to the evolution equation Eq. o is not straight¬ 
forward , since it is inte^rodifferential equation . However, as we will show below, for a special case which 
corresponds to a semiclassical approximation, this method is applicable. It is convenient to rewrite the 
wave function T in terms of the “action” S [24], 


^(y,K) = 


(2.18) 


Then the evolution equation Eq. (1.8) takes the form of a nonlinear PDE 

F{Y,l,S,j,uj) = io{Y,l) -x{ 7{Y,1)) - P {Y,l,a)) = 0 


where we introduced shorthand notations 


uj{Y,l) = 


dS [Y-l) 


7iY,l) - 1 = 


dS {Y-l) 


(2.19) 


( 2 . 20 ) 


dY ' ' ^ ' dl 

and introduced a new variable I = In (k -|- a) which remains finite in the small-K limit. The explicit form 
of the function P (Z, 7 {Y, Z, a)) in Eq. (p.20|) will be specified below in section 2.4, here we will only mention 


that P (Z, 7 {Y, Z, a)) 0. From the exact solutions of the massless BFKL I.IC , we expect that at large 


K, the effective action should depend linearly on rapidity Y and a new variable Z, 

^SC ~ ^00 T (Too 1 )^) 


( 2 . 21 ) 


where 700 and iVoo = a;BFKL( 7 oo) are constants. In a semiclassical approximation which is valid for 
the moderate values of Y, we assume that* ujsc and jsc are slowly varying functions of variables Z, Y: 


*We introduce a subscript index SC for semiclassical approximation 
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viz. u'yiYJ) « {Y,l), UJ[{Y,1) « u;2(y,/) and7^(l^,0 « {I - l {Y,l))\ i,{Y,l) « {I - ^ {Y,l)f. 

Making this assumption, the equation Eq. ( 2.19D has a form of the PDE which can be solved using the 
method of characteristics [22,23]. The characteristic lines l{t),Y(t), S{t), uit) and 7sc(i)) where t is some 
parameter (analog of time in case of classical mechanics), satisfy a system of ODEs; 


( 1 -) 

( 2 .) 


dt 07 d7sc 

dY{t) _ ^ 


dxilsc) dP{l{t),jscit),a) 


disc 


( 2 . 22 ) 


dt du) 

, , dS , , dF dF , , r 
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( 4 .) 
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dl 

dt \dl ’ dS) 

dusc (dF dF\ 

- = — ' -h W T7V7 I = 0 


di 

dP {l{t),iscit),a) 

dl 


di 


I + ^sc 


dt 


V dY 


dS J 


The second line of equation Eq. ( 2.22| ) implies that the parameter t corresponds to a rapidity Y. Erom the 
fifth line of Eq. ( p.22 ) we can see that ujsc is conserved on characteristic lines, and can be fixed from the 
asymptotic conditions Eq. (1.10|) as 


= xilsc) + Pilit),iscit),a{ioo)) = x(7oo). 


(2.23) 


A combination of Eq. (2.22)-4 and Eq. ( 2.22| )-1 allows us to eliminate the T-dependence, and find isc 
as a function of Z on a trajectory. 


disc 

dl 


dP{l{t),lsc{t),a) / (dxilsc) ^ {l{t),lsc{t),a) 


dl 


disc 


disc 


(2.24) 


Using Eq. 


-1 we can write the equation for the trajectory I = I (E), which has the form 

dl (E) dx iisc (I {Y))) dP {l{Y).lsc{l {Y)), a) 


dY 


dl 


disc 


(2.25) 


The lines I (E) = Isc (Y) give the set of trajectories. The trajectory that leads to the dominant 
contribution to ik (E, 1) can be found from the equation 


^ = lsc{Y;ioo) = lsp{Y;ioo) 


(2.26) 


Corresponding isc (Y, Isc (Y ); loo) = ISP (Y, Isp (E); i^o)- Einally, Eq. (|2.22|) -3 can be re-written in the 
form 

^ = {lsc{Y,l)-l) + x{lSc{Y,l)) + P {I, ISC (Y,I, a)) (2.27) 






































Comparing Eq. ( 2.22| ) with Eq. (|2.13 ) one can see that the effective intercept is equal to ujeff {Y,l) = 

(^s{y,Isp{y)) + izsp(y))/y. 

Before applying the semiclassical approach to the BEKL equation for massive gluons, in the following 


subsection |2.3|, we would like to test it on the massless BEKL equation (1.9), for which analytical solutions 


are known (see for example Ref. [3,6]), and only after that in subsection 2.4 we apply it to the massive 
case. 


2.3 Semi-classical approach: massless BFKL equation 


For Eq. (|L9|) E {'jsp {Y,l) ,1) = —x{l) (see Eq. (ll.lOD ). Combining this equation with Eq. (|2.12D , we 
obtain for the saddle point 75P (y, / = In/t) 


{iSP (6) + ? = 0 where ^ = l/Y. 


(2.28) 


The solution to Eq. (2.28) is shown in Fig. |3|-a. From Eq. ( ]2.13D one has that 

^BFKL (y,^ = Ihk) oc where oJeffiO = x(7sp(0)+ (0 “ (2-29) 


The dependence Wg// (0 is plotted in Fig. l-b. The Fig. j^c shows the ratio R{^,Y) defined in Eq. ( p.ldf ). 
This ratio is small at large values of Y and small thus justifying the applicability of semiclassical 
approximation in this region. 


■ysp(^) 






Fig. |-c 


Figure 3: Semi-classical solution for massless BFKL equation( see Eq. (O)): trajectories as function of ^ = l/Y 
(Fig. |a); effective Pomeron intercept (we// versus ^ (Fig. i-b); and the ratio R of Eq. (|2.14[) at fixed ^ as a function 
of y(Fig. |-c). 


2.3.1 Diffusion approximation 


As one can see from Fig. ^a at | 0 or in other words, at large y S> /, the trajectory jsP approaches a 

limiting value 75 p(oo) = Since x{l) is an analytic function near this point, it can be approximated by 


7(7) 


u)Q + D 





(2.30) 
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with Wo = 4In2 = 2.772 and D = 14('(3) = 16.822. In the approximation ( 2.30 ) the Eq. ( p.28 ) can be 
solved easily and yields 

75P = ^ - (2-31) 


2 2D 


Substituting (p.31[) into Eq. (|2.13|), we obtain 


'I'bfkl (^, Y) = (j)i. 


1 e 


1 


g 2 ^g‘^ 0 ^ 4 D’ 


2 2DJ ^ AttDY 

i.e., the semiclassical approach reproduces the diffusion approximation [ 6 ] for the BEKL equation. 

2.3.2 Double log approximation 

Eor small values of 7 ~ 0, the BEKL kernel can be approximated by 

xi-y) = 

7 


(2.32) 


In this limit, Eq. ( 2.28|) gives the trajectory 'ysp = and the wave function 


1 


=2^/^y _ 


27r ^3/2 y 


= Cpin{Y/l) 


1 


2vr/3/2 y- 1/2 


2 VFI 


(2.33) 


(2.34) 


Eq. (2.34) is the solution in the double log approximation [ 6 ]. 


2.4 Semi-classical approach: kernels of the BEKL equation with massive gluon 

Plugging the solution Eq. ( p.ll ) into Eq. (1.8), we obtain the equation for E ( 7 , 1) in the form 

= r(e'-a) + Cr(/, 7 ( 0 ,a;) - )C (1,^ (1) ,a), (2.35) 

where the emission kernel (1), a) is given by 


/C(/,7(0,«) = 


y'(l-f)= + i(l + t) + ^ A y'(i_i)2 + 2t(i + t) + (i_4»); 


:. (2.36) 


The integral over t in Eq. ( 2.36| ) can be evaluated analytically and expressed in terms of the Appel function 
Fit: 


K (17 (0,«) = F) u (t. 1 - 7 ( 0 , i t 1 -^ 

+ F (1,7 (0) Fi (7 (/), 1, i, 1 + 7 (/), t-, , 


(2.37) 


tSee Eqns. 9.180 - 9.184 in Ref. [21] 
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where 


n a I 1 — 1 

= 1 + - ± 2 
K K 


For practical reasons it is convenient to introdnce a function defined as 


= dt{P ^ - l) 
Jto 


/ dt- 
Jo 


(2.38) 


- 1 


+ 


a/( 1 - t)2 + {2/k) (1 + t) + (1 - 40)7^2 a/( 1 - t)2-l Jo 

[' dt (t-'> -1) ^ ^ ^ 

Jo 


Y^(l-t)2 + (2/K)t(l + t) + (l-4a)t2/K2 1^(1 - t)2- 


(2.39) 


Then Eq. ( 2.35 ) can be cast into the form 

E = -X(7(0) + T{l,a) + CT{l,^{l),a) - P{l,j{l),a) = -x{l{l)) - P{l, 7 il),a), (2.40) 

where 


5 1+ R - 1 1 , 

(a-1) ^B[ - ,l-7,a), 


t(/,a) = T(e^-a^ - L(e\a 
^ dt 


(2.41) 

(2.42) 


L {k, a) = 


dt 


t ^/ii^^tyTJoTt)~T^ ^Jn-tf + lt{i + t) + (i-4a)^ 

- _. _ Kin ( k(1 — k + \/—4a + (1 + k)^ j 

= — ln(l + a — K + 7/(1 — a + k)^ ) + In (l + y/l — Aa + 4k) - . - 

^ W ^ ^ ^-4a + (1 + k )2 

K In ^1 — 4a + 3 k + y^l — 4a + 4 ka/— 4a + (1 + k)^^ 


^-4a + (1 + k )2 


(2.43) 


and B {x,p,q) is the incomplete Beta function!. Fig. ^ illustrates how all the ingredients of Eq. ( 2.401) 
behave as functions of I = ln(K + a{iy)). 


3. Semi-classical solutions to the BFKL equation with massive gluon: numerical results 
3.1 Trajectories and intercepts 


For solving the general BFKL equation with massive gluons, we need to find the trajectory from Eq. ( ^.12 ) 
which will be a function of both ^ = l/Y and 1. According to the method of characteristics discussed in 
Section 2.2.2, these trajectories are the solutions of Eq. ( |2.23 ) or, equivalently, Eq. (2.24). Unfortunately, 
we are able to solve Eq. (2.23|) only numerically, and these solutions are shown in Fig. |^. 


^See Eq. 8.39 in Ref. [21] 
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1= ln(A: + a) 


7sc(0 




1 _ 



versus I = ln(/c + a). Function Ysc (0 is the inverse function to I — Isc iX) of E*!- (2-26). 


All the trajectories can be characterized by their asymptotic boundary condition 750 (0 




We note that in Eq. (|2.8f) one should understand 7 as 700 , thus reducing it to the form 

T {-a) = - X (Too) • 


Too* 


(3.1) 


m 


At small values of I all the trajectories xsc 0 vanish, reflecting the Kq (y^r)-behavior of the solution 
the coordinate representation. The negative value of xsc for soft I in a numerical solution of Eq. ( 2.231) 
is different from what one expects from a massless BEKL. In section 3.4.3 we will discuss this region in 
more detail, here we only point out that E ( 7 , 1) is an analytical function of 7 for negative values of 7 ’s 
without any singularities at 7 = —n,n = 0,1,2... . Eig. ^ shows that 75 P (Z) given by the solution of 
Eq. ( 2.23 ) and presented in Eig. |5|, satisfy the equation even at I < Igoft- 

The trajectories Isp (E) can be found by solving Eq. ( p.26 ), and are shown in Eig. |^b for several 
boundary conditions. In the ^ <C 1 region the solution may be constructed analytically without solving an 
additional equation. The first observation is that at large I, the set of trajectories should coincide with the 
same set for the massless BEKL equations (|2.28|), 


ISP {1 ,0 7™ (0 = Too- 


(3.2) 


Assuming that 75 ™^ (0 is close to ^ at small ^ (in so called Bjorken limit), we can consider a deviation 
as a small parameter. 
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wsc (0 




Figure 6: Intercept of the massive BFKL Pomeron versus I on the trajectories with different 700 • 


7SP(0 = ISC { l,7oo = -] + d'ysp{l,(). 


(3.3) 


To find 57Sp(/,0,Eq. (|^) can be re-written in the form 

(tsc (^ 7 oo = ^) + e - ^) 


hsp (0 = 1 ™ - ? 

e^O \ 


dx{'r='rsc{i,'rcx=^) + (-) dP{i,i='rsc{i,icx=^)+>^M^oa)) 

d-y dy 


(3.4) 


where e ~ 0 is a small cutoff needed to regularize some intermediate results. In particular, Eq. (3.4) 
reproduces Eq. ( 2.31 ) and gives the final result at / S> 1 where _j_ dP{i,y-ysc{i),a{y^)) ^ ^ 


dy 


The solution to Eq. ( 2.12| ) in this approximation has a form 
T(y,/) = ct>in{lSp{Y,l)) 


\ 


1 

2tt 

d^Eiysp(Y,l),l) 

Y 



e-¥ 


(3.5) 

(3.6) 


where we introduced an effective intercept 

^eff {1,0 = -E {'JSP {Y, 0,0 + ^75P (+, 0 - 2 ) 
and 75 P is given by Eq. (^^). 

The /-dependence of the effective intercepts UJ {1,0 is shown in Fig. Note at small 0 these intercepts 
are close to the massless BFKL given by Fig. ^a in the entire region of 1. In the region of large 0 the 
effective intercepts are considerably smaller than the intercept of the massless BFKL Pomeron. For small 
I, i.e., the scattering amplitude at small values of I and large values of is suppressed in comparison to 
the behavior determined by the intercept of the massless BFKL Pomeron. Since, the asymptotic region at 
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high energies (large Y) corresponds to ^ <C 1, the high energy behavior of the massive BFKL Pomeron is 
the same as for the massless one. Therefore, we have the confirmation of our basic results of the numerical 
calculation (see Ref. [1]) in the framework of the semiclassical approach. 


msp(/) 






Figure 7: uieff = uisp of the massive BFKL Pomeron versus I at fixed ^ = l/T.Fig. illustrates that ujsp (1) are 
close for small values of 


3.2 Accuracy of the semiclassical approach 


The accuracy of the semiclassical approximation is controlled by the ratio R of Eq. (2.14). Fig. ^ shows 
this ratio for different values of From the figure we conclude that for ^ < 1 the ratio R is small and we 
can safely use the semiclassical approximation. 

For I < 4, the ratio R is small in the region of large Y (small |). The reason for this is since for 
large I, the solution of the BFKL equation for the massive gluons, should coincide with the solution of the 
massless BFKL equation. These estimates confirm our expectation that the semiclassical method provides 


m i) 


R(^, 1) 




Fig. 


Figure 8: The ratio of Eq. ( ^.14 ) for large values of Y. Fig. ||-a shows this ratio at fixed ^ = l/Y which are small, 
and Fig. |b presents i? at ^ = 1 as a function of 1. 
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a reliable approach at high energies (large value of Y). 


At large our procedure does not work even at large 1. In terms of kinematics, we need to deal 
with 'ysc ^ 1, while in the Fig. ^ we used ^sc For this region of large I and we develop an 

approximation which corresponds to the double log approximation and for very large I , and coincides with 
the DLA for the massless BFKL equation. 


In Fig. ^we plot the ratio 


Rsc 


_ m _ 

{i-jsciiyf 


(3.8) 


This parameter controls the smoothness of the functions u}{Y,l) and 'y{Y,l), and thus a precision of the 
semiclassical approach. The Fig. ^ shows that this ratio is very small at large I, and doesn’t exceed ~ 0.4 
even at small 1. 


RscH) 



Figure 9: The ratio of Eq. (^) versus 1. 


3.3 Saturation momentum 


It is well known that we can find a saturation momentum by searching for the particular trajectory on 
which the wave function T (T, k) = const (see Refs. [3,25,26]). Mathematically, this corresponds to a 
solution of a system of two equations 


..ajecory: y + , = 0. 

O'J 

front line: E {'jsp {Y, 1) ,1) Y = {1 - 'ysp) I- 


(3.9) 


In the case of massless the BFKL equation, the solution to the equations of Eq. (|3.9|) is ^sp = Icr = 0.37 
[3]. For the massive BFKL equation, the critical trajectory is shown in Fig. |T^a. The equation for the 
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saturation momentum for massless BFKL equations of the form 


hr = \n{Ql{Y)/Ql{Y = ^)) = (3.10) 

1 - 7cr 

The solution to Eq. ( J3.9D for hr = In [Q^ (Y) /QliY = 0)) is shown in Fig. |^b. The difference between 
these two cases is sizable, only for small values of I = In (k + a). 


7cv(l) 



IcAY) 



Figure 10: The critical trajectories for BFKL equation with m = 0 and with m ^ 0 (see Fig. |T^a). Fig. I^b shows 
the dependence of the saturation scale (In (Q^ (T) /Ql {Y = 0))) as function of Y. 


3.4 Analytical solutions 

In this section we develop two analytical methods of searching for solutions based on the diffusion and 
DLA approximations ,for the massless BFKL equation. 


3.4.1 Diffusion approximation 


A brief glance at the trajectories for small valnes of ^ (see Fig. P) allows ns to conclude that these trajectories 
are close to 7 = | at least for I > 5. Therefore, for such values of I we can develop the diffnsion approach, 
in complete analogy with the case of massless BFKL equation, that has been discussed in section |2.3.1 . In 
the vicinity of 7 = | we can expand the general expression of Eq. ( p.40 ) as 

(3.11) 


a; ( 7 ,/) = -E{'y,l,a) = A (/) + Ai{l)[j-- 


+ A 2 (0 ( 7 - 2 


Functions A (1 ), Ai (/) and A 2 (1) are plotted in Fig. 0 -a. We see that at large I (say at I > Iq ^ 5) 
the functions Aj reach constant values, A ( 
expansion (|3.11) into Eq. (|2.25), we obtain 


the functions A* reach constant values, A (1) ^^ ■ °> wq; Ai (Z) 0; and A 2 (Z) D. Substituting 


A(Z) + Ai(Z) 


lsc{l) - 2 


+ 


A 2 (0 ( 


isp-^ 


= XI 2 


(3.12) 


16 - 
































whose solution is 

ISC (0 = 7 d (0 = ( “^1 


(0 ± ^A2(!) -4A2(i) / (2 A 2 (() 


(3.13) 


The physical corresponds to the minus sign in Eq. ( |3.13| ). In Fig. 12 we compare the trajectory Eq. ( p.l3D 
with the exact trajectory at ^ = 0.01, that has been calculated in Section 3.1, For I > Iq 10 we can 
safely use the solution Eq. (|3.13|) . 


7scil) 



Figure 12: Functions 'ysc {1^0 versus 1: exact solution to Eq. (2.23); diffusion approximation and approximation 
at small values of 7 . In the figure ^ = 0.01. 


Using Eq. (^.121) , we can calculate the wave function of Eq. (|2.13D which takes the form 

.e-¥ 


m{Y,l) = cPinhspil)) 


1 


4Z\3 (0 T 


(3.14) 
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with 


We //(^0 = xQ) - 


(? (7d(0 - I))' 


2(Ai(0 + 2A2(0 {lDil)-l))' 

The Eq. ( 3.15| ) provides a good description of the intercept for rather large Z > Zq ~ 10- 


(3.15) 


3.4.2 Small 7 approximation 

For I < Iq, we cannot use the diffusion approximation since, as we can see from Fig. [T^, the diffusion 
trajectory is much larger than the exact one in this region. Actually, the exact 7 sp(/,C) small ^ 


approaches jsp —^0. At 7 —)• 0 the general expression for P{l,j,a) in Eq. (|2.40D can be simplihed and 
takes the form 


X( 7 )+ P(/, 7 ,o) ^ -(l-e , 

7 n {l,a) 1 \ J 


1 1 


7 ' ... H{l,a)j 


with H {l,a) = 1/\J{1 + exp (—/))^ — 4a exp {—21) and /a = lna, 


The Eq. ( 2.12 ) takes the form 
1 




I 7 |p V J -fsc J 


(3.16) 


(3.17) 


The Eq. ( 3.1?| ) has two solutions in different regions: (1) 7 <C 1 but 7 (/ — la) ^ 1; and (2) 7 <C 1 
and J [I — la) <. 1- 

In the first kinematic region Eq. (S) reduces to 

1 


2 ^ ili^) 1 

Tsc 


with the solution 


ISC = \lH (I, a) 


(3.18) 


(3.19) 


We can check that 7 (/ — la) = ^ Therefore, this solution corresponds 

to the DLA approximation in this kinematic region. 

For the kinematic region 7 <C 1 and 'y {I — la) 1 the Eq. ( |3.16 ) leads to the analytical function at 
7 —^ 0 in contrast to the case of the massless BFKL kernel. Therefore, we can search for a parametrization 
w {1,0, which has the same form as Eq. ([1.12 ), viz.: 


1 ( 7,0 = -E {j,l,a = 4:) = d{l) + di{l)j + ^2 (0 7 ^, 


(3.20) 


with functions di is plotted in Fig. in^b. We can see from Fig. ^ that we can rely on Eq. ( [l.20D only for 
/-/a<l/(7 = 0.2)«5. 


- 18 - 

































Eq. (|2.12|) with cj (/,^) given by Eq. ( 3.20D , has a solution 

75c(0 =7s(0 = 1d{1) = ± -4d2(/) {d{l) /(2d2(0)- (3.21) 

The solution (p.21) is shown in Fig. 12. This solution leads to good approximation for Z = 3 -7 5 with 




^7^(0 


2(di (/) + 2d2 (0 75(0)' 


The wave function takes the form 

d/(y,o = </>*n(75p(0) 


4 da (0>" 


g-Z g<^e//(0^. 


(3.22) 


(3.23) 


3.4.3 Z < Isoft 

Earlier we have seen that the trajectory 75 p(Z) has a node at some value Isoft, and becomes even negative 
for I < Isoft- As we have seen earlier, the function utspi'y) is analytic at 7 = 0 , so we can extrapolate 
to negative but small 7 ’s using Eq. (^.2Clf ) (see Fig. |l^. In principle, we can expect some high twists 
singularities at 7 = —n with n = 0 , 1 , 2 ,... which stem from the expansion of 


y^(l - t)2 + (2 /k) (1 +t) + (1 - Aa)/^? 


Y, Cn{l,a) f‘ 


(3.24) 


71=0 


at small t in Eq. (2.39). Taking the integral over t in Eq. ( p.39 ) in the vicinity of small t, one can see that 
for 7 ^ —n we obtain the same expression as in Eq. ( p.l 6 | ), replacing 7 in Eq. ( p.l 6 ) by 7 + n. The energy 
(intercept) turns out to be a regular function at 7 = —n, and we can use Eq. (3.20|) for I 7 + 1| < 1 with 


function {1) calculated at 7 = —1 ( see Fig. 12 and Fig. 13). 


Fig. 1^ illustrates that the intercept ui ( 7 , 1) is the analytical function without singularities at negative 
7 which increases at large | 7 |. 

Unfortunately, we have not found a simple analytical approach that is able to describe the scattering 


amplitude at all values of 1. However, we would like to recall that the numerical solutions to Eq. (3.3) and 
Eq. ( |3.4| ) depend neither on the value of the QCD coupling, nor on the initial condition, and reduce the 
procedure of calculation of the scattering amplitude to a simple equation. Solving this equation is much 
easier task than the exact numerical calculation of the eigenvalues and eigenfunctions that was done in 
Ref. [1]. 


4. Conclusions 

In our previous paper (see Ref. [1]) we studied the BFKL equation with massive gluons in the lattice and 
proved that its spectrum coincides with the spectrum of the massless BFKL equation. This observation 
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Figure 13: Function w ( 7 , 1) versus 7 at fixed 1. lobfkl = X (^)- 


gives rise to a hope that the correct large impact parameter (6) behavior of the scattering amplitude 
A oc exp(—m6), which is the inherent feature of the massive BFKL equation, will not affect the high 
energy behavior of the scattering amplitude. Therefore, we can expect that the modification of the BFKL 
equation due to confinement would not affect strongly the main equations that govern the physics at high 
energy ( in particular, the non-linear equations of the CGC/saturation approach to high density QCD). 


In this paper we developed the semiclassical approximation which allowed us to investigate the high 
energy behavior of the scattering amplitude. The method provides a simple procedure for the calculation 
and reduces it to a numerical solution of Eq. ( p.23 ), which is much simpler than the direct numerical 
calculations of the eigenvalue problem in the lattice realized in Ref. [1]. 


Having these solutions, we propose a modification of high energy asymptotic behavior, caused by the 
correct large b exponential fall off of the amplitude. Actually, we did not find any unexpected behavior, 
and the semiclassical solution reproduces the scattering amplitude which is very close to the amplitude of 
the massless BFKL equation, at least at high energies. 

In section |3.3| we estimate the value of the saturation momentum, solving the linear evolution equation 
with very general assumptions about the non-linear corrections. We demonstrated that the value of the 
saturation momentum is close to the one for the massless BFKL equations, leading to the assumption that 
saturation physics will look similar for both massive and massless BFKL equations. 


We believe that in this paper we taken the natural next step in the understanding of the influence of 
the correct large b decrease of the amplitude on its high energy behavior. It should be stressed that this 
behavior is interesting as it provides the high energy amplitude for the electroweak-weak theory, which 
can be measured experimentally. The solution which has been discussed in this paper, determines the 
asymptotic high energy behavior of the electroweak-weak theory for zero Weinberg angle. We plan to 
address the physical case of nonzero Weinberg angle [2] elsewhere. 
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